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Preface 

This  paper  describes  a  solution  to  the  problem  of 
optimally  controlling  the  drift  of  a  space  vehicle  stationed 
at  one  of  the  earth-moon  triangular  llbratlon  points.  Many 
1 1  bra 1 1  on  point  studies  have  preceded  this  one  and  I  am 
greatly  Indebted  to  the  foundation  which  they  have  provided. 
Of  particular  Importance  to  this  study  were  the  Investiga¬ 
tions  by  Dr.  Lynn  E.  Wolaver  and  other  studies  performed 
under  Dr.  Wolaver's  sponsorship,  as  described  In  Chapter  I 
of  this  report. 

This  thesis  began  as  a  follow-on  study  to  the  analog 
computer  Investigation  by  Captain  Gerald  T.  Rudolph  (Ref  2). 
However,  after  much  effort,  It  became  evident  that  using  a 
trial  and  error  approach  In  the  style  of  Captain  Rudolph, 

I  could  not  hope  to  obtain  satisfactory  results.  This 
approach  was  finally  abandoned  In  favor  of  modern  optimal 
control  theory  lechnlques.  The  change  was  made  with  reser¬ 
vation  due  to  the  belief  that  a  satisfactory  formulation  and 
solution  In  this  manner  might  not  he  possible.  Fortunately, 

I  was  able  to  formulate  the  problem  In  a  manner  which  made 
It  readily  solvable  on  a  high-speed  digital  computer. 

I  wish  to  express  my  appreciation  to  my  thesis  advisor. 
Captain  Edward  A.  Kern,  for  his  support  and  advice.  I  also 
want  to  thank  Dr.  Lynn  E.  Ifolaver  of  the  Aerospace  Research 
Laboratories,  who  proposed  and  sponsored  this  thesis,  for 
his  help  and  suggestions. 
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Abstract 


A  control  system  Is  devised  for  maintaining  a  space 
vehicle  In  close  proximity  to  one  of  the  earth-moon  tri¬ 
angular  llbratlon  points  while  minimizing  fuel  consumption. 
The  problem  Is  formulated  as  an  optimal  state  regulator 
problem  of  variational  calculus  and  then  modern  control 
theory  Is  applied.  A  quadratic  performance  criterion  Is 
used  which  leads  to  a  linear  feedback  control  system.  The 
feedback  gains  are  obtained  by  solving  a  matrix  differential 
equation  of  the  Rleattl  type  on  a  high-speed  digital 
computer.  Performance  of  the  resulting  optimal  system  Is 
verified  and  further  analyzed  on  an  analog  computer. 
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OPTIMAL  STATION  KEEPING  AT 
THE  L4  L I Q RATION  POINT 

I .  I  introduction 

Because  of  the  strong  perturbing  Influence  of  the  sun, 
a  space  vehicle  stationed  at  one  of  the  triangular  Lagrange 
llbratlon  points  of  the  earth-moon  system  would  require  a 
thrust  control  system  to  confine  Its  drift.  It  Is  the 
object  of  this  thesis  to  devise  a  realizable  control  system 
which  uses  a  minimum  of  fuel  to  restrict  the  drift  of  a 
space  vehicle  stationed  near  one  of  these  llbratlon  points. 

Background 

The  Very  Restricted  Four-Body  Problem.  The  llbratlon 
points  are  equilibrium  points  at  which  a  p;  tide  would,  if 
given  the  proper  Initial  velocity,  remain  fixed  relative  to 
the  two  other  masses  of  a  restricted  three-body  system.  The 
adjective  "restricted"  Imposes  two  conditions;  the  particle's 
mass  must  be  so  small  that  It  does  not  affect  the  motion  of 
the  two  larger  masses  and  the  two  larger  masses  must  move 
in  circular  orbits  about  their  common  center  of  mass. 
Lagrange,  in  1772  proved  that  five  such  points  exist.  Their 
locations,  fixed  in  a  rotating  coordinate  system,  are  illus¬ 
trated  in  Figure  1.  It  is  well  established  for  this  three- 
body  system  that  if  the  mass  ratio  of  the  two  large  bodies 
is  less  than  0.04,  then  L4  and  L5  are  stable  In  the  sense 
that  a  small  disturbance  results  in  a  non-di  vergent  motion 
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Figure  1 

L I  BRAT  I  ON  POINTS  OF  THE  RESTRICTED  THREE-BODY  PROBLEM 
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close  to  the  llbratlon  point.  The  other  three  equilibrium 
points  are  unstable  for  all  mass  ratios. 

For  the  case  of  the  earth-moon  system,  the  effect  of 
a  fourth  body,  the  sun,  can  not  be  ignored.  A  mathematical 
model  which  has  been  used  successfully  to  represent  the 
earth-moon-sun  system  is  the  model  known  as  "the  very 
restricted  four-body  problem"  (Refs  2,  3,  4,  &  5).  This 
model  retains  the  conditions  of  the  restricted  three-body 
problem  and  Imposes  the  additional  conditions  that  the 
center  of  mass  of  the  earth-moon  system  move  In  a  circular 
orbit  about  the  center  of  mass  of  the  entire  system  and  fhat 
the  orbits  be  coplaner  (see  Figure  2).  Equations  of  motion 
for  this  model  have  been  derived  and  nondlmensl onal 1  zed 
(Refs  3  and  4)  by  setting  the  unit  of  distance  equal  to  the 
earth- to-moon  distance,  the  unit  of  mass  equal  to  the  total 
mass  of  the  earth  and  moon,  and  the  unit  of  time  such  that 
uj  Is  unity  and  one  time  unit  Is  equal  to  4.348  mean  solar 
days.  The  nondlmenslonal  equations  are 

x-2y-x-(Ms/R2)cos$=3P/3x  (1) 

y  +  2 x  -  y  +  ( M s / R 2 )  sin  $  *  3P/3y  (2) 


in  which 


m 


1  -  m 


P  = 


[(x  -  1  +  m)2  +  y2y*  [(x  +  m)2  +  y2]'2 


C ( x  +  !  cos  ip)2  +  (y  -  R  sin  $)2I! 
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Figure  2 

VERY  RESTRICTED  FOUR-BODY  REPRESENTATION 
OF  THE  EARTH-MOON-SUN  SYSTEM 
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$  ■  wt  +  a 

m  ■  Mm  (In  nondl mens  1 ona 1  mass  units) 
w  *  (u>  -  fi )  /w 

When  linearized  about  the  old  L4  llbratlon  point  the  equation 
becomes  (Refs  5:80-81  and  4:37-39) 

x  -  2 y  -  CjX  -  C2y  B  C 4  +  C6  cos  2$  +  C7  sin  2* 

+  C8  cos  «j>  +  C g  sin  <p 
+  Cj  jn  (x  cos  2 4>  -  y  sin  2* ) 

+  C12n  ( 3 X l  cos  <p  -  Yj_  sin  <j>)x 

+  C !  2  n  (Yl  cos  <j>  -  sin  4> )  y  (3) 

••  • 

y  +  2x  -  C2x  -  C3y  »  C5  +  C7  cos  2$  -  C6  sin  2<j> 

-  Cg  cos  <j>  +  C1Q  sin  <t> 

-  Cj  j  n  (x  sin  2<p  +  y  cos  2 4> ) 

+  C12n  (YL  cos  4>  -  X l  sin  <p )  x 

+  C12n  ( X |_  cos  <f>  -  3Yl  sin  <t>)y  (4) 

in  which  (in  normalized  units) 

XL  =  (1  -  2m)/2  C7  =  -  3nYL/2 

Yl  =/f/2  C8  «  -  3 n ( 3 X l  +  Y£)/8R 

n  =  Ms/R3  Cg  =  3nXLYL/4R 

Cj  =  (3  +  2n)/4  C10  =  3n(3Y^  +  X[)/8R 

C2  =3/3(1-  2m)/4  CL1  =  3/2 

C3  =  (9  +  2n)/4  C12  =  -  3/4 R 

C4  *  nXL/2 

C5  =  nYL/2 

C6  *  3nXL/2  (5) 
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Minimizing  Drift  by  Selection  o  f  I  nit  1  a  1  Conditions.  I  n 
a  Master?  Thesis  by  Capt.  Paul  M.  Ulshafer  (Ref  4)  and  in  a 
study  by  Dr.  Lynn  F.  Wolaver  (Ref  5),  Initial  condition?  were 
derived  which  minimized  the  drift  of  a  space  vehicle 
stationed  near  the  L4  point.  The  studies  were  based  on  the 
model  given  by  the  very  restricted  four-body  problem.  It 
was  found  that  by  making  a  judicious  choice  of  initial 
conditions,  which  include  the  initial  sun  direction  a,  the 
drift  from  the  L4  point  remains  less  than  0.013  distance 
units  during  a  one-year  mission. 

Ana  1 oq  Simulation.  In  a  Masters  Thesis  by  Capt.  Gerald 
T.  Rudolph  (Ref  2),  an  attempt  was  made  using  an  analog 
computer  to  "minimize  the  energy  required  to  hold  a  space 
vehicle  at  the  earth-moon  L4  libration  point."  The  study 
failed  to  accomplish  this  stated  Intention,  however,  it  did 
demonstrate  the  following: 

(1)  the  feasibility  of  using  an  analog  computer  to 
simulate  the  equations  of  motion. 

(2)  the  accuracy  of  using  the  linearized  equations,  Eqs 

(3)  and  (4),  without  the  terms  involving 

C8,  Cgt  Cjq,  and  C12. 

Known  Cost.  Prior  to  this  investigation,  only  the  two 
extremes  of  L4  station  keeping  cost  were  known,  that  for  no 
thrust  control  and  that  for  perfect  control.  Without  thrust 
control,  thrust  cost  is  of  course  zero;  however,  the  cost  in 
drift  would  be,  assuming  optimum  initial  conditions,  0.013 
for  a  one-year  mission.  Perfect  control,  i.e.,  continuous 
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maintenance  of  the  satellite  precisely  at  L4,  would  have  no 
cost  in  drift,  however,  the  cost  In  thrust  would  be  finite 
for  any  given  mission  duration.  The  thrust  acceleration 
required  for  perfect  control  would  be  the  amount  necessary 
to  exactly  cancel  the  forcing  function  of  Eqs  (3)  and  (4), 

1  .e. , 


ux  -  -  (C,  +  C6  cos  2 4>  +  C7  sin  2<f> 

+  C8  cos  $  +  C9  s  1  n  $ )  (6) 

Uy  -  -  (Cs  +  C7  cos  2$  -  C6  sin  2$ 

-C9cos<|>  +  C10sinif>)  (7) 

The  average  magnitude  of  the  required  thrust  acceleration 
would  be 

uav  '  <>/2»)  +  uy>%  d*  (8) 

To  compute  Uay  analytically  would  be  difficult,  therefore 
numeral  integration  using  the  digital  computer  was  employed. 
The  terms  containing  C0 ,  C9,  and  C10  were  neglected  giving 
an  approximate  result  (in  normalized  units)  of 

Uay  =  0.008576  (9) 

Between  the  two  extremes  of  permitting  no  drift  and 
permitting  an  uncontrolled  drift,  the  control  cost  was  unknown 
prior  to  this  study.  For  all  that  was  known  at  the  time,  the 
optimum  thrust  policy  might  have  been  to  maintain  perfect 
contro 1 . 
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Problem  Formu I ati on 

Perf ormance  I ndex .  The  purpose  of  this  thesis  is  to 
minimize  the  fuel  requirements  of  a  space  vehicle  stationed 
near  the  earth-moon  L4  llbration  point.  Although  only  the 
L4  point  is  considered,  the  results  apply  equally  well  to 
the  15  point  because  of  symmetry. 

In  order  that  the  study  results  be  independent  of 
vehicle  weight  or  rocket  performance,  total  velocity  Incre¬ 
ment,  AV,  Is  used  as  a  cost  measure  rather  than  fuel. 

AV  is  calculated  from  the  relation 

t 

A V  -  /  ||  U(t )  ||  dt  (10) 

o 

in  which  jJ(t)  is  the  acceleration  due  to  the  thrust.  aV 
Is  related  to  fuel  according  to  the  well-known  rocket 
equation 

Fuel  Weight  »  [Final  Vehicle  Weight]  [exp  ( a V/ g0 1 sp) -  1  ]  ( 1 1 ) 

Therefore,  if  it  is  assumed  that  Isp  is  constant  and  that 
no  staging  occurs,  minimum  AV  implies  minimum  fuel.  The 
performance  index  is  then  the  aV  required  to  hold  the 
satellite  within  some  given  peak  drift  distance  for  a  given 
amount  of  time.  No  consideration  will  be  given  to  attitude 
control  costs. 

System  Equations.  The  model  of  the  Very  Restricted 
Four-Body  Problem  is  assumed  in  this  study.  If  the  terms 
involving  C8,  C9,  C10»  and  Cj2  are  dropped  from  the 
linearized  equations  and  a  thrust  acceleration  term  is  added. 
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nondlmenslonallzed  system  equations  of  sufficient  accuracy 

are  .jbtalned; 

i 

I 

i  x  -  2y  -  C  j  x  -  C2y  ■  C4  +  C6  cos  2$  +  C7  sin  2  <(. 

+  cun  (x  cos  2<p  -  y  sin  2$)  +  Ux  (12) 

y  +  2x  -  C2x  -  C3y  ■  C5  +  C?  cos  2$  ••  C6  sin  2<*> 

-C11n(xsin2<ti+ycos2<|))  +  Uy  (13) 

The  coefficients  were  evaluated  based  on  the  astronomical 

'  constants  (in  normalized  units) 

1/m  «  82.301 

MS 

n 

ve suiting  in  the  values 

w  =  0.9252 
XL  =  0.4378 
Yl  =  0.3660 
n  =  0.005595 
Ci  =  0.7528 
C2  *  1.267 
C3  «  2.253 
C4  =  0.001365 
C  5  =  0.002423 

C6  =  0.004094 
C7  =  -  0.007268 
Cu  =  3/2 


=  328,905.2 

»  0.07480133  (14) 
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Unless  otherwise  noted,  all  quantities  are  expressed 
in  normalized  units,  The  normalized  units  can  be  converted 
into  physical  units  using  the  following  relations: 

1  time  unit  -  4.348  mean  solar  days 
1  distance  unit  *  23C>,0Q0  statute  miles 
1  velocity  unit  =  3360  ft/sec 

1  acceleration  unit  =  0.008943  ft/sec2  (16) 


Approach 

It  was  believed  at  first  that  a  trial  and  error  approach 
using  an  analog  computer  would  yield  a  solution  to  the 
problem.  Several  thrusting  schemes  were  postulated  and 
tested  which  did  actually  limit  the  drift  to  any  desired 
value  for  any  duration.  However,  measurement  of  the  cost 
(aV)  indicated  no  saving  when  compared  to  the  cost  of 
perfect  control.  One  of  the  schemes  te^ed  used  a  control 
design  based  on  a  root  locus  analysis  oi  .ne  transfer 
function  x/Uy.  Beautiful  stability  resulted,  however,  cost 
measurements  could  not  prove  its  merit. 

The  problem  was  successfully  solved  by  considering  it 
as  an  optimal  regulator  problem  of  variational  calculus.  The 
direct  application  of  modern  control  theory  is  at  best  a 
difficult  task  due  to  several  peculiarities  of  the  problem. 

These  are: 

(1)  the  constraint  of  peak  rift, 

(2)  the  time  varing  coefficients  of  the  system 
equations  and 

(3)  the  absence  of  a  terminal  cost  or  end-point  conditions. 
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The  approach  used,  therefore,  was  to  first  devise  a  substitute 
performance  criterion  which  could  be  handled  by  Pontryagl n 1 s 
maximum  principle,  solve  for  the  optimal  control,  and  finally 
measure  its  worth  In  terms  of  the  actual  performance  criterion. 
It  was  decided  to  use,  without  constraints,  the  quadratic  cost 
functi on 

J  -  l  J1  [q  (x2  +  y2)  +  U2  +  U2]dt  (17) 

1  o  y 

hopeful  that  perhaps  the  drift  could  be  controlled  to  any 
desired  peak  value  merely  by  varying  the  weighing  factor  q. 
This  particular  performance  index  was  chosen  for  two  very 
Important  reasons: 

(1)  It  leads  to  a  readily  solvable  control  law 

(2)  the  control  turns  out  to  be  a  linear  function 

of  state,  and  hence  a  closed-loop  control  system 
resul ts . 

A  closed-loop  control  system  is  highly  desirable  since  the 
system  equations  do  not  take  into  account  the  real-world 
perturbations  or  even  the  nonlinearities  of  the  very 
restricted  four-body  problem.  With  a  feedback  (closed  loop) 
control  it  becomes  more  likely  that  such  perturbations  have 
little  effect. 

The  remainder  of  this  report  describes  the  technique 
and  results  of  taking  this  final  approach  to  the  problem. 
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I  I .  Analysis 


Theory 

For  a  linear  system,  the  optimal  control  with  respect 
to  a  quadratic  criterion  can  be  obtained  from  a  rather 
eloquent  theory.  It  Is  eloquent  not  only  from  a  mathematical 
point  of  view  but  also  from  an  engineering  point  of  view,  for 
the  optimal  control  turns  out  to  be  a  linear  function  of 
state.  The  theory  was  developed  by  Kalman,  but  a  derivation 
based  on  Pontryagin's  maximum  principle  can  be  found  in 
Optimal  Control  by  Athans  and  Falb  (Ref  1:750-766).  A 
summary  of  the  control  law  which  results  from  this  theory 
follows  (Ref  1 : 762-763) ; 


Given  the  linear  system 

X(t)  =  A(t)X(t)  +  B(t)U(t)  (18) 

and  the  cost  functional 
0  a  2  <  X.(  T )  ,  F  X(T)> 

+  5  /  [<X(t),  Q(t)X(t)>  +  <U(t),  R(t)U(t)>]dt  (19) 

2  t0 

where  U_(t)  is  not  constrained,  T  is  specified, 

F  and  Q(t)  are  positive  semi  definite,  and  R(t) 
is  positive  definite.  Then  an  optimal  control 
exists,  is  unique,  and  is  given  by  the  equation 


U(t)  =  -  R_1(t)B'(t)K(t)X(t)  (20) 


where  the  n  x  n  symmetric  matrix  K(t)  is  the 
unique  solution  of  the  Riccati  equation 
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K(t)  -  -  K ( t )  A  ( t )  ••  A  '  ( t )  K ( t ) 

+  K(t)B(t)R”  1  (  t)B‘  (  L)K(t)  -  Q(  t.)  (?1) 

satisfying  the  boundary  condition 

K(T)  -  F  (2?) 

The  state  of  the  optimal  system  is  then  the 
solution  of  the  linear  differential  equation 

X(t)  -  CA(t)  -  B(t)R"1 (t)B' (t)K(t)1  X(t)  (23) 

Figure  3  shows  the  structure  of  the  optimal  system. 

Note  that  the  optimal  control  is  generated  by  a  linear,  time- 
varying,  feedback  gain.  Since  R(t)  nd  B(t)  are  known 
matrices.  It  follows  that  the  matrix  Kt  )  determines  the 
gain  of  the  feedback  path  and  the  behavior  of  the  system. 

The  matrix  K(t)  is  therefore  referred  to  as  the  "gain" 
matrix. 

The  gain  matrix  is  the  solution  of  a  nonlinear  matrix 
equation,  Eq  (21),  and  for  this  reason  a  closed  form  solution 
is  not  usually  obtainable.  This  is  no  great  obstacle  how¬ 
ever,  since  a  numerical  solution  of  the  gain  matrix  can  be 
obtained  easily  using  a  digital  computer.  Note  that  Eq  (21) 
must  be  solved  backward  in  time  from  time  T  in  order  that 
the  boundary  condition  given  by  Eq  (22)  is  satisfied. 

One  fortunate  result  of  the  theory  is  that  the  gain 
matrix  K(t)  is  independent  of  the  state.  Once  the  system, 
the  cost  J,  and  the  terminal  time  T  have  been  specified, 
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K(t)  can  be  computed.  Initial  conditions  of  the  states 
have  no  effect  on  the  design  of  the  optimal  control  system. 
The  Initial  conditions  affect  only  the  cost. 


Application 

The  system  equations  used  in  this  analysis  were  presented 
In  the  introduction  (see  page  9  ).  "or  ease  of  reference  they 
are  repeated  here  in  a  slightly  different  form: 


x 


2y  +  ax  +  by  +  h  (x  cos  2$  -  y  sin  2$) 


+  d  +  f  cos  2 4>  -  g  sin  2<j>  +  Ux  (24) 

y  3  -  2x  +  bx  +  cy  -  h  (x  sin  2<s>  +  y  $  2$) 


where 


+  e  -  g  cos  2<j>  -  f  sin  2<j>  +  (25) 

0  3  wt  +  a 

a  3  initial  direction  of  the  sun 
w  =  0.9252 
a  3  Cj  3  0.7528 
b  3  C2  =  1  .267 
c  =  C3  =  2.253 
d  =  CM  =  0.001365 
e  =  C5  =  0.002423 
f  3  C6  =  T. 004094 
g  =  -C7  ==  0.007268 

h  =  Cj  jn  3  0.008393  (  26) 
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Before  the  control  low  can  be  applied,  the  system 
equations  must  be  put  In  the  form  of  Eq  (1°),  that  Is  of 
the  form 

X(t)  ■  A(t)X(t)  +  B(t)U(t)  (18) 

It  Is  not  clear,  at  first  glance,  that  this  task  Is  possible 
due  to  the  time  varying  forcing  functions  which  appear  In 
the  system  equations.  However,  by  defining  an  extra  three 
state  variables  (X5,  X6,  and  X7),  the  necessary  formulation 
can  be  accompl 1  shed . 

The  state  variables  are  defined  as  follows 

Xt  »  x 
X2  »  x 
X  3  a  y 
xu  -  y 

x5  •  1 

X 6  =  cos  2$ 

X7  ■  sin  2<p  (27) 

They  are  Inter-related  by  the  following  equations: 

ii  =  x  2 

X2  =  2X4  +  (a  +  h  cos  2<t>)Xj 

+  (b  -  h  sin  2 4> )  X 3  +  dX5 
+  fX6  -  gX7  +  Ux 

X3  ■  Xu 


1  6 


t 

i 
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■  -  2  X2  +  (b  -  h  sin  2$ ) X l 

+  (c  -  h  cos  20  )X3  +  oX5 
-  gX6  -  fX7  +  Uy 

X5  ■  0 
X6  «  -  2wX7 

X,  =  2wX6  (28) 


The  system  equations  are  now  expressable  In  the  required 
form,  i.e. 


in  which 


X(t)  =  A(t)X(t)  +  B(t)U(t) 


X(t) 


Xj(t) 

x2(t) 

x3(t) 

x4(t) 

x5(t) 

x6(t) 

X7(t) 


(18) 


(29) 


U(t) 


ux(t) 

Uy(t) 


(30) 
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A  ( t ) 


n 

1 

0 

0 

0 

0 

0 

(a  +  h  cos 

2*) 

0 

(t>  - 

h  sin 

2^) 

2 

d 

f 

-g 

0 

0 

0 

1 

0 

0 

0 

( b  ■  h  sin 

2  $ ) 

-2 

(c  - 

h  cos 

20) 

0 

e 

-g 

-  i* 

0 

0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

0 

0 

-  2w 

0 

0 

0 

0 

0 

2  w 

0 

(31  ) 


B(t)  =  B 


0  0 
1  0 
0  0 
0  1 
0  0 
0  0 
0  0 


(32) 


The  cost  functional  used  in  this  analysis,  as  discussed 
previously  (see  page  11)  is  the  quadratic  cost 


J  =  \  /T  Cq  (X2  +  Xn  +  (U 2  +  U2)]dt  (33) 

<-  o  *  y 

which  can  be  expressed  in  the  form  of  Eq  (19) 

J  =  \  <X(T),  F  X(T)>  +  \  fT  C<x(t),  q( t)x( t)> 

*o 

+  <U(t)  ,  R(t)U(  t)>]dt  (19) 
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by  letting 


F  ■  COD  a  7  by  7  zero  matrix 


(34) 


Q(t) 


q 

o 

o 

o 

o 

0 

0 


0 

0 

0 

0 

0 

0 

0 


0 

0 

q 

o 

o 

o 

o 


o 

o 

o 

o 

o 

0 

0 


0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 


(35) 


R(t)  =  R  * 


1 

0 


0 

1 


(36) 


t0  =  0 


(37) 


Equations  (20)  and  (23)  can  now  be  simplied  somewhat  by 
carrying  out  the  indicated  matrix  multiplications.  The 
following  equations  result: 


U(t)  =  - 


'  2  1 


'4  1 


'2  2 


'  4  2 


'2  3 


'4  3 


'2  4 


>4  4 


'2  5 


'4  5 


'2  6 


'4  6 


'2  7 


•4  7 


X(t) 

-1  (38) 
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/ 


— 

- 

\ 

0 

0 

0 

0 

0 

0 

0 

K2  i 

K?  2 

K2  3 

K2  4 

K  2  r 

K2  6 

K2  7 

0 

0 

0 

0 

0 

0 

0 

<  Aft)  - 

Km 

K4  2 

K4  3 

K44 

K4  5 

K4  6 

K4  7 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 

- 

y 

)  x(t) 


(39) 


where  the  gain  matrix  is  defined  as 


K(t)  = 


Ki  i  K i  2 

•/  y 

*  2  1  Ts  2  2 


K7 1  K7 2 


•  •  • 


K 1  7 

Kz  7 


•  •  • 


K7  7 


(40) 


This  completes  the  adaptation  of  the  problem  to  fit  the 
control  law  stated  in  the  previous  section  (see  page  12). 

Note  that  the  problem  now  satisfies  all  the  necessary  conditions 
required  by  the  control  law.  Eq  (21)  must  now  be  solved  (by 
some  numerical  technique)  for  the  elements  of  the  gain  matrix. 
Then  Eqs  (38)  and  (39)  can  be  used  to  obtain  the  optimal 
control  and  the  system  performance. 
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III.  Digital  Computati ons 

The  nature  of  the  gain  matrix  and  the  performance  of  the 
system  is  investigated  in  this  section.  The  gain  matrix  and 
the  system  performance  was  computed  with  the  aic  of  an  IBM 
7094  digital  computer  using  the  Runge-Kutta  method  of 
numerical  Integration. 

In  the  course  of  this  investigation  it  becomes  evident 
that  modifications  can  be  made  in  the  application  of  the 
control  law.  These  modifications  are  motivated  by  the  desire 
to  reduce  the  complexity  of  the  control  system.  Fortunately, 
several  other  considerations  lend  support  to  the  same  modi¬ 
fications,  namely  considerations  of  the  ease  of  computation, 
the  ease  of  making  an  analog  simulation,  and  improved  per¬ 
formance  in  terms  of  the  actual  cost  of  the  system  (see 
page  8  ). 

The  Gal n  Matrix 

Figure  4  is  a  sketch  which  illustrates  the  nature  of  the 
gain  matrix  elements.  In  order  to  show  all  characteristics 
in  a  single  illustration,  a  composite  of  the  gain  matrix 
elements  is  depicted  rather  than  any  one  particular  element. 
Each  gain  matrix  element  displays  the  same  general  trend, 
that  of  remaining  at  some  steady  state  level  with  a  slight 
oscillation  until  the  terminal  time  T  is  approached,  then 
decaying  to  zero.  The  elements  differ  in  the  oscillation 
amplitude  in  the  steady  state  region,  the  point  at  which  the 
decay  begins,  and  the  character  (smooth  or  oscillatory)  of 
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the  decay. 

Three  traits  of  the  gain  elements  become  loss  pro¬ 
nounced  as  q,  weighting  on  drift,  Is  Increased.  These 
traits  are:  oscillations  during  the  steady  state  interval, 
peaks  during  the  transient  interval,  and  a  relatively  long 
duration  transient  Interval.  As  q  Is  increased,  each  of 
the  gain  elements  exhibit  more  nearly  a  stop  function 
character,  as  viewed  backwards  In  time.  That  is,  the 
oscillation  amplitudes  decrease,  the  peaks  in  the  transient 
Interval  diminish,  and  the  transient  interval  becomes  very 
small.  For  a  q  of  0.1,  the  most  severe  oscillation  is 
only  about  6%  (ak  relative  to  k),  exhibited  by  element  K2  6 , 
This  element  also  exhibits  the  most  severe  transient  interval, 
about  twenty  normalized  time  units  for  a  q  of  0.1.  Since 
the  oscillations  are,  in  general,  of  such  small  magnitudes 
and  because  the  transient  intervals  are  generally  short 
relative  to  a  one-year  mission  (84  normalized  time  units),  a 
modification  of  the  control  law  is  considered  which  ignores 
the  transient  interval  and  the  oscillation.  The  modified 
control  law  differs  from  the  optimal  law  only  in  that  it 
uses  a  constant  gain  matrix  rather  than  the  time-varying  one. 
The  constant  gain  matrix  is  set  equal  to  the  average  steady- 
state  level  of  the  time-varying  gain  matrix.  System  per¬ 
formance  using  inis  modified  control  law  is  described  in  the 
next  section. 

The  constant  gain  matrix  is  somewhat  less  difficult  to 
compute  than  the  time-varying  gain  matrix.  Also,  the  compu¬ 
tation  no  longer  depends  upon  the  terminal  time  T.  A 
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constant  "A"  matrix  whose  elements  are  the  average  values  of 
the  time-varying  A  matrix  elements  Is  used  for  the  computation. 
The  time-varying  elements  of  the  A  matrix  are 


A2  j 

■  0.7528 

+  0.008393 

cos 

2* 

a2  3 

-  1.267 

-  0.008393 

s  1  n 

2* 

A*,  \ 

-  A2  3 

-  2.253 

-  0.008393 

cos 

2  4> 

(41) 

The  approximation  Is  therefore  a  very  close  one  since  the 
time-varying  values  never  differ  from  their  average  values 
by  more  than  0.008393.  The  results  of  calculating  the  gain 
matrix  In  this  fashion  are  tabulated  in  Tables  I  and  II  for 
values  of  q  ranging  from  .01  to  1,000,000.  Only  those 
elements  required  by  Eqs  (38)  and  (39)  are  listed.  The 
computer  program  used  to  compute  this  data  Is  described  In 
Appendix  A. 

System  Performance 

System  performance  using  the  modified  control  law  was 
Investigated  and  compared  with  performance  of  the  optimal 
system.  Results  of  this  comparison  show  that  the  modified 
control  system  is  superior  to  the  optimal  system  when  actual 
costs  are  considered  rather  than  quadratic  cost. 

System  performance  data  was  computed  for  the  optimal 
control  and  the  modified  control  systems  using  a  q  of  1.0 
and  a  mission  time,  of  15  (about  two  months).  This  data  is 
listed,  for  the  last  0.4  time  units  of  the  mission,  in 
Tables  III  and  IV.  Peak  drift  with  the  optimal  control 
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Table  II 

STEADY-STATE  GAINS  (CONTINUATION) 
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Table  III 

PERFORMANCE  OF  THE  OPTIMAL  CONTROL 

INITIAL 

CONDI TIONS 

XI  *  0. 

X2»  0. 

X3"  U. 

X  4  =  0. 

SUN  DIRECTION* 

0.  DEGREES 

Time 

Control 

Drift 

Peak  Drift 

aV  Cost 

14.6000 

4.5  7906E-Q4 

4. 0524  IE-03 

4. 83045E-03 

6.48640E-02 

14.6100 

4. 36649E-04 

4.09946E-03 

4.83045E-0  3 

6.48684E-02 

14.6200 

4. 15838E-04 

4.  14693E-03 

4.83045E-03 

6.48725E-02 

14.6300 

3.9  5479E-04 

4.  194856,-03 

4.8  3045E-03 

6, 48765E -02 

14.6400 

3. 75574E-04 

4. 24326E-03 

4.83045E-03 

6.48802E-02 

14.6500 

3 .561285-04 

4.2S221E-03 

4.83045E-03 

6.48838E-02 

14.6600 

3. 371466-04 

4. 341 72E-03 

4.8  3045E-03 

6. 488726-02 

14.6700 

3.  186316-04 

4. 391 84E-03 

4.830456-03 

6. 48904E -02 

14.6600 

3.00589E-04 

4.44261E-03 

4.83045E-03 

6.48934E -02 

14.6900 

2.8  3022E-04 

4.49407E -03 

4.83045E-03 

6.48962E-0? 

14.7000 

2.659356-04 

4. 54627E-03 

4.8  3045E-0  3 

6.48989E-02 

14.7100 

2.49332E-04 

4. 59924E-03 

4.83045E-03 

6.4901 3E -02 

14. 7200 

2. 3321 7 E -04 

4.65304E-03 

4.83045E-0  3 

6.4903  7E -02 

14.7300 

2  •  17595E -04 

4.70770E-03 

4.83045E-03 

6.49059E -02 

14.7400 

2. 02470E-04 

4. 763266-03 

4.830456-03 

6.49079E-02 

14.7500 

1 . 87845E-04 

4. 81 977E-03 

4.83045E-03 

6.49098E-02 

14,7600 

1.73726E-04 

4.87727E-03 

4.87727E-03 

6.4911 5E-02 

14.7700 

1.6011 7E- 04 

4.935806-03 

4.93580E-03 

6.49131E -02 

14.7800 

1.47023E-04 

4,  995406-03 

4.99540E-03 

6.49146E -02 

14.7900 

1 . 34447E -04 

5.05612E-03 

5.05  6 12E-0  3 

6.491596-02 

14.8000 

1  •  2239  56-04 

5. 117996-03 

5.1 1799E-03 

6.4917  IE -02 

14.8100 

l  ,  10871 E -04 

5. 181056-03 

5  •  1 8 105E-0  3 

6.49182E -02 

14.8200 

9.98815E-05 

5.245346-03 

5.24534E-03 

6.49192E-02 

14.8300 

8. 9430  2  E-0  5 

5 . 3  1090E-0  3 

5.31090E-0  3 

6. 49201E -02 

14.8400 

7.952296-05 

5. 37777E -03 

5.37777E-03 

6.492096-02 

14.8500 

7. 01648E-05 

5. 445976-03 

5.4459  7F- 03 

6. 492166-02 

14.8600 

6.13617E-05 

5. 51 555E-03 

5.51555E-03 

6.49223E-02 

14.8700 

5.311936-05 

5.586546-03 

5  .5  8654E-03 

6.492286-02 

14.8800 

4, 54436E-05 

5.65896E-03 

5.658966-03 

6.49232E-02 

14.8900 

3. 8340  76-05 

5. 732856-03 

5.7  328  5E-03 

6 • 49236E -02 

14.9000 

3. 18170E-0  5 

5.8O824E-03 

5.8O824E-03 

6.49239E-02 

14.9100 

2. 587916-05 

5.8851 56-03 

5.88515E-03 

6  .49242E-02 

14.9200 

2.053386-U5 

5 . 96  360E-0  3 

5.96  360E-03 

6.492446-02 

14.9300 

1.57882E-05 

6.04362E-03 

6.04362E-03 

6.49246E-02 

14.9400 

1.  16494E-05 

6.  12524E-03 

6.12524E-0  3 

6. 492476-02 

14.9500 

8. 12499E-06 

6.208466-03 

6 .208466-03 

6. 492486-02 

14.9600 

5.2227  86-06 

6.  29  33  IE- 0  3 

6.29331E-03 

6. 49248E-02 

14.9700 

2  •  95076E -0 6 

6.  37 98  IE -0 3 

6.379816-03 

6.492486-02 

14.9800 

1.317226-06 

6. 46796E-03 

6.46  7 9 6 1 - 0 3 

6.492496-02 

14.9900 

3.307316-07 

<  557786-03 

6.55  77  8E-03 

6.492496-02 

15.0000 

0. 

6. 649296-03 

6.649296-03 

6.49249E-02 
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Table  IV 

PERFORMANCE  OF  THE  MODI 

FILD  CONTROL 

INITIAL  CONDITIONS 

XI-  O. 

X3 *  0 . 

SUN  DIRECT  ION» 

X  0  . 

X4=  0. 

0.  DEGREES 

T 1  mo 

Con  tro 1 

Drift 

Peak  Drift 

aV  Cost 

14  .6000 

4.30926E-03 

3.  7802  3E-03 

4.839286-03 

6. 770605 -02 

14.6130 

4 • 34  13  16-03 

3. 814306-03 

4.83928E-03 

6.78294E-02 

14.6200 

4.37  29  96 -03 

3.  848106-03 

4.B3928E-03 

6.78731E-02 

14  .6300 

4.  404266-03 

3. 881586-03 

4.839286-03 

* • 79171c -02 

14.6430 

4.4351  IE-03 

3.914736-03 

4.839286-03 

6.796156*02 

14.6530 

4. 465496-03 

3.94752F-03 

4. 839286-03 

6. 800616-02 

14.6600 

4.  49  539F-03 

3.  97993E-03 

4.839286-03 

6.80511E-02 

14  .6700 

4.  52 47 7E -03 

4.01 192F-03 

4.83928E-03 

6.80963E-02 

14.6800 

4. 55  36  0E-03 

4.04  349F -03 

4.839286-03 

6.81419E-02 

14.6900 

4. 5818  76-03 

4.074606-03 

4.839286-03 

6.818775-02 

14.7000 

4. 609546-03 

4.  105246-03 

4.83923E-03 

6.82338E-02 

14.7100 

4.53  66  OE-O  3 

4.  135396-03 

4. 83928E-03 

6.82802E-02 

14.7200 

4.663026-03 

4. 165036-03 

4.83928E-03 

6.  83268E -02 

14.7300 

4. 6887  8E-03 

4.  19414E-03 

4.83928E-03 

6.83737E-02 

14.7400 

4.  713866-03 

4.  222706-03 

4.839286-03 

6.842086 -02 

14.7500 

4.738246-03 

4. 250706-03 

4. 83928E-03 

6.846826-02 

14.7600 

4.761906-03 

4.27813E-03 

4.83928E-03 

6.851  58E-02 

14.7700 

4. 784836-03 

4. 304966-03 

4.83928E-03 

6.856376-02 

14.7800 

4.  807006-03 

4. 331196-03 

4.83928E-03 

6.861175 -02 

14.7900 

4. 828406-03 

4. 35680E-03 

4.83928E-03 

6. 86600E -02 

14.8000 

4.849016-03 

4.38177E-03 

4.B3928E-03 

6.870855-02 

14.8100 

4. 868836-03 

4.406116-03 

4.83928E-03 

6.875725-02 

14.8230 

4.887836-03 

4. 42979E-03 

4.83928E-03 

6. 88061E-02 

14.8330 

4. 906016-03 

4.452816-03 

4.83928E-03 

6.885516-02 

14.8400 

4.  92335E-03 

4. 47516E-03 

4.8  3928E-03 

6. 89044E -02 

14.8500 

4. 939846-03 

4.  496826-03 

4.83928E-03 

6.895385-02 

14.8600 

4. 95548E-03 

4.517806-03 

4.H3928E-03 

6. 90033E -02 

14.8730 

4.9702 4E-03 

4. 53808E-03 

4.83920E-O3 

6.905306-02 

14.8800 

4.984!.  36-03 

4. 55766E-03 

4 .3  3928E-03 

6.910295-02 

14.8900 

4.99713E-03 

4. 57653E-03 

4.83928E-03 

6.915285-02 

14.9330 

5.009246-03 

4. 59469E-03 

4.83928E-03 

6. 92029E -02 

14.9100 

5. 02045E-03 

4.612136-03 

4.83928E-03 

6.92531E -02 

14  .9200 

5.030750-03 

4. 628846-03 

4. 339286-03 

6.930346-02 

14.9330 

5.04015E-03 

4.644846-03 

4. 83928E-03 

6. 935386-02 

14.9400 

5. 04863E-03 

4.660106-03 

4.R3928E-03 

6. 94043E -02 

14.9500 

5.3562DE-Q3 

4.674636-03 

4.83928E-03 

6.945495 -02 

14  .9630 

5.  J6285E-03 

4. 68844E-03 

4.33928E-03 

6.950555-02 

14.9700 

5.  06857E-03 

4.70  15  IE-03 

4.3  3928  E-0  3 

6.95562E-02 

14.9800 

5.07338E-03 

4.  71 38  5E -  03 

4.83928E-03 

6. 96069E -02 

14.9900 

5.07726E-03 

4.725456-03 

4.83928E-03 

6. 965775-02 

15.0330 

5.08022E-03 

4. 736336-03 

4.83928E-03 

6.97085E -02 
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remains  constant  at  0.0  0-183  until  the  last  0,25  time  units 
of  the  mission,  then  Increases  to  a  value  of  O.OOG65  at  the 
terminal  time  of  15.  With  the  modified  control,  peak  drift 
never  exceeds  0.00484  during  the  entire  mission.  Cost  (aV) 
Is  about  7%  greater  (0.0697  verses  0.0649)  than  ttiat  of  the 
optimal  control  system  but  peak  drift  is  smaller  by  more 
than  25%  (0.00484  verses  0.00665).  For  longer  missions, 
the  difference  in  cost  would  be  even  less,  while  the 
reduction  In  peak  drift  would  be  exactly  the  same.  The 
comparison  Is  summarized  in  Table  V. 


Table  V 

COMPARISON  -  OPTIMAL  VS  MODIFIED 

CONTROL  ( T-l 5) 

Optimal  Control 

Modified  Concrol 

/  (drift  squared)  dt 

0.0002153 

0.0002304 

/  (control  squared)  dt 

0.0003088 

0.0003370 

Total  Quadratic  Cost 

0.0005241 

0.0005674 

Peak  Drift 

0.006649 

0 . 004 S3  9 

AV  Cost 

0.06492 

0.06971 

Performance  calculations  on  the  modified  control  system 
were  made  using  a  number  of  q  values  which  range  from  0.01 
to  10.  The  computer  program  used  is  described  in  Appendix  C. 

As  a  result  of  two  important  findings,  the  performance 
of  the  modified  control  system  can  be  summarized  in  an 
astonishingly  simple  graph  which  is  shown  in  Figure  b.  The 
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findings  which  permit  this  type  of  display  are  that 

(1)  After  a  br<ef  transitory  period,  the  peak  drift  is 
uniquely  determined  by  q  (neither  mission  du¬ 
ration  nor  initial  conditions  affect  the  drift 
after  the  transient) . 

(2)  After  the  transitory  period,  the  control  becomes 
periodic  with  u  magnitude  which  is  roughly  constant 
(the  mean  value  of  this  control  magnitude  therefore 
can  be  used  to  represent  the  system  cost  in  terms 
of  cost  per  unit  time). 

The  graph  of  Figure  5  essentially  answers  the  question 
posed  earlier  (see  page  7  )  of  whether  or  not  it  would  be 
less  costly  to  thrust,  keeping  the  satellite  on  station 
perfectly  rather  than  permitting  a  finite  drift.  As  seen 
by  the  graph,  cost  decreases,  almost  linearly,  with  permissi¬ 
ble  drift.  Perfect  station  keeping  is  decidedly  a  more 
costly  option. 

Further  study  of  the  modified  control  system  was 
accomplished  with  the  aid  of  an  analog  computer.  This  ad¬ 
ditional  study  is  described  in  the  next  chapter  along  with 
more  details  of  the  oerformance  traits  noted  in  this  chapter. 
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I  V  .  At i n  1  ruj  S  i  mu  I  <i  H  on 

The  system  dynamics  and  tin?  modified  control  system  were 
simulated  on  an  Electronic  Associates,  Inc.  "PACT  TR48"  analog 
computer.  This  machine  has  a  repetitive  operation  feature  and 
a  built-in  oscilloscope  which,  when  used  together,  provides  a 
near  instantaneous  display  of  the  system  response.  Recordings 
were  made  on  an  "X-Y"  plotter  of  the  most  illustrative  and 
significant  runs  (with  the  computer  in  normal  operation  mode). 

The  control  is  generated  by  1A  constant  gain  feedback 
paths.  Two  control  systems  were  wired  and  tested,  one  having 
gains  appropriate  for  a  q  of  0.01  and  the  other  for  a  q 

of  1.0  (see  Tables  I  and  II  for  1 1.  gain  values).  The 
circuitry  used  for  the  simulation  is  described  in  Appendix  C. 

Several  important  aspects  of  the  system  are  discovered 
which  would  have  been  difficult  to  detect  from  digital  data. 
Discoveries  were  made  concerning  the  steady-state  tendency, 
the  effect  of  initial  conditions,  and  the  existence  of  a 
stable  orbit  without  control. 

Steady-State  Tendency 

Regardless  of  initial  conditions  or  the  value  of  q, 
the  control  system  always  maneuvers  the  satellite  into  a 
cyclic,  steady-state  trajectory  about  the  LA  point.  The 
period  of  rotation  about  the  LA  point  in  this  state  is  always 
precisely  half  that  of  the  sun  (it/w  time  units).  This 
steady-state  tendency  is  illustrated  in  Figure  6.  The 
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steady- sta to  t.  raj  or  tory  resembles  an  off-center  ellipse. 

I  *  r,  size  is  a  function  only  of  the  value  of  ci  .  For  compari¬ 
son  purposes,  the  trajectory  for  a  q  of  zero  (no  control) 
is  illustrated  in  Figure  7.  Once  the  satellite  reaches  the 
steady-state  path,  it  remains  there  for  as  long  as  the  control 
system  is  in  operation. 

Although  Initial  conditions  have  no  effect  on  the  steady- 
state  trajectory,  they  do,  however,  affect  the  initial  path 
of  the  satellite.  lie  fore  reaching  the  steady-state  conc'ition, 
depending  upon  initial  conditions,  the  satellite  may  drift 
out  even  beyond  the  peak  drift  of  the  steady-state  trajectory 
as  in  Finure  6.  (luring  the  transient,  the  control  system  is 
apparently  manuevering  the  satellite  to  place  it  on  the 
steady-state  path  with  the  proper  velocity  and  a  particular 
phase  relationship  relative  to  the  time  of  month.  Initial 
conditions  can  always  be  chosen,  of  course,  which  would  place 
the  satellite  precisely  on  the  steady-state  trajectory 
initially,  therefore  eliminating  the  transient  path  entirely. 
For  a  given  time  of  month,  i.e.,  initial  sun  direction  a, 
such  a  set  of  initial  conditions  are  obviously  unique. 

Initial  Conditions  for  Mi n i m u m  Cost 

Once  the  steady -state  conditions  are  established,  the 
control  becomes  cyclic  much  the  same  as  does  the  satellite 
trajectory.  The  control  is  plotted  as  a  vector  ( Uy  verses 
Ux  with  time  as  a  parameter)  in  Figure  R  for  a  q  of  1.0 
arid  with  zero  initial  conditions.  During  the  steady-state 
hase,  the  control  magnitude  is  roughly  constant  and  t hero- 
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Figure  7 

TRAJECTORY  WITHOUT  CONTROL 


35 


J 


GRE/GE/67- TO 


C.C.C/EE/G7  10 


fore,  tho  system  cost  Increases  linearly  with  time  during 
this  phase.  The  magnitude  of  tho  control  during  the  steady- 
state  phase  Is  unaffected  b\  Initial  conditions,  however 
this  Is  not  so  during  the  transient  phase.  Tho  magnitude  of 
the  control  during  the  transient  phase,  and  therefore  tho 
system  cost,  can  be  minimized  by  a  proper  choice  of  Initial 
conditions.  Such  a  set  and  Its  effect  on  the  control  Is 
Illustrated  In  Figure  9.  Although  tho  control  can  be  mini¬ 
mized,  It  Is  done  so  at  the  sacrifice  of  permitting  the 
satellite  to  drift  beyond  the  steady-state  path.  A  compro¬ 
mise  Is  possible  however,  as  shown  In  Figure  10.  In  this 
situation  the  satellite  remains  completely  Inside  the  steady- 
state  path  but  the  reduction  In  the  control  magnitude  Is  not 
quite  as  great  as  before.  For  a  q  of  0.01,  no  better 
Initial  conditions  are  found  than  those  which  place  the 
satellite  on  the  steady-state  trajectory  initially  (see 
Figures  11  and  12). 

No  matter  how  much  the  control  is  minimized  during  the 
transient  phase,  the  effect  on  the  system  cost  will  be  small 
if  long  duration  missions  are  considered.  In  Figure  13 
the  system  cost  with  a  minimized  control  magnitude  during  the 
transient  period  is  compared  with  the  system  cost  in  which 
there  is  no  transient  period.  Certainly  the  saving  is 
insignificant  for  missions  of  more  than  a  few  months. 

Optimum  Initial  Conditions  Wi thout  Control 

Extrapolation  of  the  "cost  verses  peak  drift"  curve  of 
Fi  ire  5  suggests  that  a  stable  orbit  inving  a  peak  drift  of 
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Figure  9 

RESPONSE  WITH  "MINIMUM  COST 
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Initial  Conditions: 
a  ■  129.2  degrees 
x  ■  -0.00061 
x  »  -0,01442 
y  -  -0.00524 
y  ■  0.00444 
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about  0.011  would  exist  If  q  were  zero,  that  Is  to  say 
without  a  control  system  at  all!  Further  Investigation 
shows  that,  for  any  value  of  q,  the  controlled  system 
under  steady-state  conditions  always  exhibits  the  same  phase 
relationships  between  sun  direction,  satellite  position  and 
velocity  (see  Figures  14  thru  17).  This  fact  Implies  that 
perhaps  a  particular  phase  relationship  Is  mandatory  for  the 
establishment  of  stable  orbits,  and  In  particular  for  a 

stable  orbit  with  a  q  of  zero  (no  control).  Indeed  this 

turns  out  to  be  the  case.  If  the  same  phase  relations  are 
established  by  appropriate  selection  of  the  Initial  con¬ 
ditions,  the  satellite  does  In  fact  remain  In  a  stable  orbit 
without  a  control  system  and  furthermore  one  that  has  a  peak 
drift  of  about  0,011.  This  Is  not  Inconsistent  with  the 

findings  of  Dr.  Lynn  E.  Wolver  (Ref  5)  and  Capt.  Paul 

Ulshafer  (Ref  4)  which  show  that  a  fairly  stable  orbit  can 
be  achieved  having  a  peak  drift  on  only  0.013  during  a 
one-year  mission  provided  a  judicious  choice  of  Initial  con¬ 
ditions  is  made. 

The  appropriate  initial  conditions  can  be  derived  from 
data  of  systems  with  small  q  values.  The  required  data 
are  the  values  of  Xj ,  X2,  X3,  X4,  and  X7  as  functions  of 
time  after  reaching  steady-state  conditions.  A  time  plot 
of  these  variables  in  steady-state  conditions  is  shown  in 
Figures  18  and  19  for  a  q  of  0.01.  The  optimum  initial 
conditions  for  any  given  a  are  those  which  correspond  to 
an  X7  having  the  value  sin  ?.a  since 
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X7  ■  sin  ?.  (w  t  +  a)  (27) 

Fair  accuracy  would  be  obtained  using  the  data  for  a  system 
with  a  q  of  0.01,  but  for  a  more  accurate  determination, 
the  data  from  Tables  VI  (q  *  0.01),  VII  (q  ■  0.02)  and 
VIII  (q  ■  0.03)  can  be  extrapolated  to  give  values  which 
would  correspond  to  a  q  of  zero. 

,  v 

The  analog  plot  of  Figure  20  Illustrates  the  stable 
orbit  without  control.  For  this  plot,  the  extrapolation 
method  just  described  was  used  to  determine  the  Initial 
conditions.  The  Initial  sun  direction  ct  was  chosen  such 
that  X7  corresponded  to  an  exact  table  entry  (7.2707599  E - 0 3 ) 
but  any  a  could  have  been  chosen.  The  orbit  remains 
reasonably  stable  for  nearly  three  years  according  to  the 
analog  plot.  The  slight  instability  Is  due  only  to  accuracy 
limitations  of  the  analog  computer.  The  potentiometer 
settings  become* extremel y  critical  for  uncontrolled  orbits  of 
long  duratl on.  such  as  the  one  depicted  here. 
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Initial  Conditions 
■  0.418  dog ree  s 
«  -0.00642' 

-  0.01391 
=  0.00685 

-  0.00496 


First  Year 
0  <  t  <  84 


Second  Year 
84  <  t  <  168 


Figure  20 

"NO-CONTROL"  TRAJECTORY 
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Model  Deficiencies 

The  model  used  In  this  study  Is  a  gross  simplification 
of  the  real  world.  The  model  does  not  account  for  such 
things  as  eccentricity  of  the  moon's  orbit,  the  pull  of 
gravity  from  other  planets,  solar  wind,  etc.  Even  the 
nonl Inearl  ties  of  the  very  restricted  four-body  problem 
were  removed  by  a  linearization  of  the  equations.  However, 
since  the  control  Is  generated  from  feedback  of  the  state 
variables,  It  Is  felt  that  these  model  deficiencies  are  not 
serl ous .  > 

Computational  Approximations 

....  T  '  .  .  ..  , 

.Approximations  were  made  on  the '.linearized  equations 

of  motion  which  aided  the  computational  process.  All 
coeff  Iclen t s  were  rounded.- off  to  four  slgnl  f  1  can t  f  1  g ures 
and  terms  Involving  coefficients  less  than  0.00002  were 

i . •  *•  ■  '  ‘*fr  V  ■  ■  ■- .  ...  .  ■  A. 

neglected.  The  theory  did  not  require  that  such  approximations 

be  made.  State  variables  for  cos  <p  and  sin  <t>  could  have 

(■ 

been  defined  which  would  handle  the  neglected  terms.  How¬ 
ever,  the  previous  work  of  Rudolph  (Ref  2)  proved  that  the 
approximations  cause  only  minute  error,  undectable  on  an 
analog  simulation. 

Real  1 zabll 1 ty  of  the  Control  System 

Thrust  Imp! Imentat i on .  The  thrust  levels  required  for 
a  typical  light-weight  satellite  are  extremely  small.  As  an 
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example,  a  1  000-lb.  satellite  requires  an  average  thrust  of 
about  0,0018  lbs.  to  confine  Its  drift  to  a  50n  mile  radius 
(see  Figure  5).  Ion  engines  or  monopropellant  rockets  would 
be  likely  candidate  for  Impllmentlng  such  low  thrust  levels. 
Higher  thrust  engines  could  also  be  used,  If  fired  briefly 
at  discrete  time  Intervals.  The  thrust  vector  rotates  with 
a  period  of  about  15  days,  therefore  only  two  firings  a  day 
might  approximate  the  continuous  thrust  vector  with  sufficient 
accuracy. 

Rudolph  proposed  (Ref  2 : 49 , 74 )  that  the  necessary  thrust 
be  obtained  from  solar  wind  pressure.  This  scheme  could  not 
succeed,  however,  because  the  solar  wind  rotates  (with 
respect  to  the  ear*  -moon  coordinate  system)  at  a  frequency 
which  Is  half  that  of  the  required  thrust  vector. 

Sensor  Requirements.  The  control  system  requires 

continuous  knowledge  of  all  seven  state  variables.  This  is 

not  considered  to  be  a  difficult  task  because  four  of  the 
•  • 

states  (x,  y,  x,  and  y)  merely  represent  position  and 
velocity.  Two  others  (sin  241  and  cos  241)  can  be  de¬ 
termined  by  sensing  the  sun's  direction.  The  only  other 
state  (Xs)  is  the  constant,  unity,  however,  the  feedback 
gains  associated  with  X5,  as  well  as  all  other  feedback 
gains,  must  be  calculated  a  priori,  based  on  the  best  known 
values  of  astronomical  constants. 

A 1 1 1 tude  Control .  Although  attitude  control  costs  were 
not  considered  In  this  study,  It  Is  conceivable  that  a 
significant  amount  of  fuel  could  be  used  to  control  attitude. 
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If  the  satellite  mission  does  not  have  Its  own  attitude 
requirements,  It  might  be  possible  to  approximate  the  proper 
rocket-nozzle  orientation  by  establishing  a  spin  of  4t 
radians  per  month,  which  is  the  rotation  froquency  of  the 
thrust  vector.  If  this  scheme  Is  used,  control  actlvatlo. 
would  be  required  only  occasionally,  If  at  all.  The  control 
would  be  needed  merely  to  establish  the  Initial  spin  and 
maintain  synchronization. 
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VI .  Summary  and  Recommendations 

ummary 

A  position  control  system  was  devised  for  a  space 
vehicle  stationed  at  the  earth-moon  L4  1 1  bra 1 1  on  point.  The 
procedure  was  to  assume  the  linearized  equations  of  the  very 
restricted  four-body  problem  and  then  apply  optimal  control 
theory  to  a  quadratic  criterion  The  optimal  control  was 
then  modified,  resulting  In  a  control  system  which  not  only 
was  less  complex  but  also  one  that  performed  better  with 
respect  to  the  more  meaningful  criterion  aV. 

in  studying  the  behavior  of  the  controlled  satellite, 
several  Interesting  discoveries  were  made: 

(1)  The  path  about  the  L4  point  becomes  periodic 
after  i  brief  transient. 

(2)  Initial  conditions  affect  performance  only 
during  the  transient. 

(3)  A  stable  orbit  can  be  obtained  without  control 
by  choosing  proper  Initial  conditions. 

Recommendations 

Further  study  Is  recommended  to  test  the  control  system 
using  a  more  realistic  model  of  the  real  world.  One  possi¬ 
bility  is  to  use  the  complete  nonlinear  equations  of  the  very 
restricted  four-body  problem.  Another  possibility,  and  one 
wi  h  would  be  a  far  more  complete  test,  is  to  use  actual 
ephemeris  data  of  the  solar  system. 

Several  other,  less  important,  studies  are  recommended: 
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(1)  Determine  the  effect  of  using  Impulse 
thrusting  to  approximate  continuous  thrust. 

(2)  Determine  the  effect  of  er" oneous  or  noisy 
sensor  Information. 

(3)  Determine  the  effect  of  using  Inaccurate 
values  of  the  astronomical  constants. 

(4)  Determine  the  effect  of  solar  wind. 
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Appendix  A 

Fortran  Program  for  Gain  Cal culatlons 

Listed  here  Is  the  digital  program,  written. In  Fortran 
IV,  which  solves  for  the  gain  matrix  K.  The  program  assumes 
a  constant  A  matrix  and  outputs  only  the  steady-state  value 
of  the  K  matrix.  The  gain  matrix  Is  computed  for  ten  differ¬ 
ent  values  of  q  beginning  with  any  specified  value  end 
Incremented  by  any  specified  amount. 

A  time  Interval  MG  must  be  specified  a  priori,  suf¬ 
ficiently  large  to  Insure  that  the  steady-state  value  of  K 
Is  reached.  For  five-place  accuracy  (five  significant 
figures),  MG  must  be  at  least  50  for  a  q  of  0.01,  25  for 
a  q  of  1.0,  and  12  for  a  q  of  10.  The  value  of  the  K 
matrix  one  time  unit  prior  to  the  end  time  Is  outputed  as  a 
check  to  Insure  the  steady-state  value  has  been  reached. 

The  14  elements  of  the  K  matrix  required  for  performance 
calculation;  are  outputed  on  punch  cards  In  the  format 
required  by  the  programs  of  Appendix  B.  The  required  Input 
data  Is: 

AA,  AB,  AC,  AD,  AE,  AF,  AG,  AT,  OM 

the  constants  a,  b,  c,  d,  e,  f,  g,  h,  and 
w  respectively 

Q  f  i  rs  t  q  va  1  ue 

DQ  q  Increment 

MG  time  Interval 

The  format  of  the  input  data  must  be  as  Illustrated  at  the  end 
of  the  program. 
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Gain  Matrix  Calculation 


*JOP  0.19.750  67-040.STEINBERG.AFI T  SE 

$  I B  JOB 

•6IRFTC  MAIN  M94.XR7 

COMMON  /COMl/A(7.7) »0,GG( 7.7) .AT.OM.ALF 
DIMENSION  Z ( 28 ) »X { 7  ) 

EXTERNAL  FA 

NAMELIST  /Nl/AA.AB. AC.AD.AE. AF.AG.AT.OM 
NAMELIST  /N4/0.DQ.MG 
READ  <  5  .  N 1  ) 

READ  (5.NA) 

DT*-1. 

DO  6  1*1,7 
DO  6  J*1 ,7 
6  A ( I > J 1 *0, 

A  (  1 , 2  )  *  1  • 

A ( 2 ♦ 1 ) ■ AA 
A ( 2 » 3 ) *AB 
A  (  2  » A  )  *  2  • 

A (2 , 5 ) 8 AD 
A ( 2  » 6 ) *  AF 
A  (  2  »  7  ) 8  —  AG 
A  (  3 .4 ) "  1  « 

A ( 4 , 1 ) 8 AB 
A<4,2)«-2. 

A ( 4  » 3 ) 8 AC 
A  (  4 , 5  )  8  A  E 
A ( 4 ♦ 6 ) =- AG 
A(4.7)=-AF 
A(6»7)»-2.*OM 
A(7,6)=2.»OM 
DO  52  NO8 1.10 
T»0. 

DO  12  I  8  1  »  28 

12  zm-o. 

NT  =  0 
14  K  8  1 

DO  16  L  8 1 ♦ 4 
DO  16  J  =  L , 7 
GG(L.J)=Z(K) 

GG< J,L)=Z(K) 

16  K=K+1 

IF  (NT.LT.MG-1 )  GO  TO  71 

WRITE  (6,13)  Q» (GG( 2. J ) , J=1 ,7 ) , ( GG( 4, J ) , J81 ,7 ) 

13  FOPMAT  (F10.4/(E17.8) ) 

IF  (NT.EO.MG)  GO  TO  72 

71  CAU.  RKDES  ( T.Z.28.DT.FA) 

NT=NT4 1 

GO  TO  14 

72  PUNCH  13,  Q»(GG(2»J)»J81,7)»(GG(4»J)»Js1»7) 
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QaQi-DO 

S2  CONTI  NUF 
STOP 
END 

SIRFTC  GAIN  M94*XR7 

SUBROUTINE  FA  (  T  » 7.  *  D  Z  ) 

COMMON  / COM  1 / A ( 7  »  7 ) »Q»  GG ( 7  *  7 ) »AT»OM»ALF 

DIMENSION  Z<  23  )  ,DZ1 28) *ZZ( 7,7 ) *ZA( 7,7) *ATZ<  7.7) »ZBZ< 7. 7 ) 

<'-l 

DC  7  1=1*7 
DO  7  J= I ,7 
7.ZC  I  ♦  J)  =  Z(K) 

zztJ»n«z(Ki 

7  K  =  K  +  1 

DO  3  1=1*7 
DO  7  Ja  I  ,7 
Z  A  (  I  »  J  >  '■»  0  • 

ATZ ( I ,J) =0. 

ZBZ( I *  J  )  =  ZZ ( I *2)*ZZ<2»J)+ZZ( I *4)*ZZ(4*J) 

DO  3  K= 1  *  7 

Z  A  (  I »  J ) =ZZ ( I  *  K ) *A  f  K  *  J ) +Z  A (  I  *  J ) 

3  ATZ ( I  » J ) =A ( K  » I  >*ZZ( K* J)+ATZ( I  * J ) 

K»1 

DO  11  1*1 *7 
DO  11  J»I  .7 

DZ(K)=-ZA<  I  *  J ) ~ATZ (  I  *  J ) +ZBZ (  I*J) 

11  K=K+ 1 

DZ(  l)=DZm-Q 
DZ I  14 ) =DZ ( 14 ) -Q 
RETURN 
END 

HIBFTC  RUNKUT  M94*XR7 

SUBROUTINE  RKDES ( X ♦ Y . N *DX ♦ F ) 

DIMENSION  Y< 28 ) ,Y0< 28 ) *YT ( 28 ) »YP( 28 ) ,P0( 28 ) 

DIMENSION  PI ( 28 ) *P2 ( 28) »P3 ( 28  ) 

1  xo  =  x 

X=X+DX 

H=DX 

2  I F ( ARS(H) .GT.ABS(X-XO) )  H=X-XO 
DO  4  1  =  1  ,N 

4  YO ( I  ) = Y (  I  ) 
ht  =  h 

XT  =  XO 

R  M  A  X  D  =  1.E37 
DO  5  1  =  1  ,N 

5  YT (  I  ) =  YO (  I  ) 

ASSIGN  6  TO  K 
GO  TO  20 
DO  7  1  =  1  ,N 
YPII >=Y(  I  ) 

KT=0.5*H 
ASSIGN  9  TO  K 
GO  TO  2  G 
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o  DO  10  I*hN 
in  YTIII>Yil) 

XT-XO+HT 
ASSIGN  11  TO  < 

20  CALL  F(XT.YT.PO) 

DO  21  r ■ 1  *  N 

21  Y(  I  )»YT(I)+0.5*HT*P0m 
CALL  F( XT+0.3*HT,Y»P1 ) 

DO  22  I - 1 *N 

22  YU)"YT(  n+,5#HT*Pl  (1) 

CALL  F(XT+0.5*HT,Y,P2) 

DO  23  I  ■  1  *N 

23  Y( I)«YT( J)+HT#P2(  I) 

CALL  F ( XT+HT  » Y  *P3  ) 

DO  24  I  » 1  »N 

24  Y(  I  )*YTm+HT<MPO(  I  )+2.*(Pl(  I  )+P2(  I  )  )+P3(  I  )  )/6. 

GO  TO  K, (6, 9,11) 

11  RMAX-0. 

DO  12  I»1,N 

12  RMAX-AMAX1 (RMAX,«07*ABS( <Y< I) -YP C I ) ) /Y ( I ) > ) 

IF  < (RMAX.GT.l. E-06), AND. IRMAX.LT. RMAXP))  GO  TO  17 
XO-XO+H 

IF ( XO.EQ.X )  RETURN 

IF ( IRHAX.LT. l.E-7). OR. < RMAX.GT. RMAXP ) )  H*H+H 
GO  TO  2 

17  H«HT  •; 

XT-XO 

DO  19  I  *1  *N 

18  YP ( I ) *YT (  I  ) 

19  YT(I)»Y0(I) 

RMAXP  «  RMAX 
GO  TO  8 

END 

SDATA  * 

tNl  AA-. 7528, AB»1. 267, AC-2, 253* AD-,00 1365, 

AE». 002423,  AF». 004094,  AG». 007268,  AT*. 008393, 
OM».92S2  $ 

SN4  Q« 1 . *  DO“ 1 •  ,  MG“2  5  $ 

$EOF 
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Appendix  D 

Fortran  Program  for  Performance  Calculatl ons 

Two  programs  are  listed  here  which  calculate  performance, 
They  differ  only  In  the  type  of  data  which  is  outputed,  The 
first  program  outputs  cost  data,  l.e.,  aV,  peak  drift,  con¬ 
trol  and  drift.  The  second  program  outputs  the  states 
X7,  X j ,  X2 ,  X3  and  X4 .  Both  programs  calculate  performance 
for  a  one-year  mission.  The  two  programs  are  basically  the 
same  and  could  easily  be  combined  Into  a  single  program. 

The  required  Input  data  for  either  program  Is: 

AA,  AB,  AC,  AD,  AE,  AF,  AG,  AT,  OM 

the  constants  a,  b,  c,  d,  e,  f,  g,  h,  and 
w  respectively 

X ( 1 ) ,  X(2),  X(3) ,  X(4) 

the  initial  values  of  states  X3,  X2 ,  X3  and 
X4  respectively 

ALF  the  Initial  sun  direction 

LM  the  number  of  q‘s  for  which  the  performance 
Is  to  be  calculated  (LM  Is  restricted  to  ten 
or  less) 

LM  sets  of  gains  as  outputed  from  the  program 
of  Appendix  A 

The  format  of  the  Input  data  must  be  as  Illustrated  at  the 
ends  of  the  two  programs. 
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Performance  Calculation  -  Costs 


SJOR  0,15,5000  67-114, STEINBERG, AFIt'sE 

SIR JOB 

$IBFTC  MAIN  M94,XR7 

COMMON  /C0M1/A(7,7) ,  GG ( 7 ♦ 7 ) , AT »0M , ALF 
DIMENSION  Z(28).X(7),Q(10)»G3(4,7,10)»XX(4) 
DIMENSION  V2( 841 ) » RP2 ( 841 ) *R2 ( 84  1 ) , T2 ( 34 1 ) ,UM( 841) 
DOUBLE  PRECISION  DV 
EXTERNAL  FX 

NAMELIST  /N1/AA,AB,AC,AD,AE,AF,AG»AT,0M 
NAMELIST  /N3/X  ,  ALF  »LM 
READ  t  5 ,N  1  ) 

READ  (5,N3) 

DO  99  J ° 1  , 4 
99  XX ( J ) - X ( J ) 

LP-LM+1 
0 ( LP  )  a0  • 

DO  43  K-2*4,2 
DO  43  J» 1 » 7 
43  G3(K,J,LP)"0« 

DO  42  NQ»1,LM 

READ  (5,13)  Q(NQ)  ,  (  (G3(K,J,NQ)  ,J-1,'7>  ,<-2,4, 2) 

42  WR I TE ( 6  » 1 3 )  0( NO ) , ( (G3(K,J,NQ) * J» 1 »  7 ) »Ka2 ,4 *  2  ) 

WR  I  TE  (  6  ♦  13  )  Q  (  LP  )  ,  (  (  G3  (  K ,  J  ,  LP  )  »  J®  1 , 7  ) .» K“2  » 4 , 2  ) 

13  FORMAT,(F10.4/(E17.8> ) 

ALFD-ALF 

ALF»ALF*3. 1415927/180* 

DO  6  I - 1  *  7 
DO  6  J-1,7 
6  A ( I , J ) *0 • 

A  (  1 , 2  )  °  1  • 

A  ,  2 » 1 )  *  AA  •  •  .  . 

A  (  2 , 3 ) °  AB 

A ( 2 ,4 ) =2  •  ; 

A ( 2 » 5 ) *AD 
A ( 2 , 6 ) * AF 
A ( 2 , 7 ) “-AG 
A ( 3 ,4 ) » 1 » 

A ( 4 , 1 ) » AB 
A ( 4 » 2 ) a-2 • 

A  (  4 , 3  ) 55  AC 
.  A ( 4 , 5 ) * AE 
A ( 4 , 6 ) a-AG 
A ( 4 , 7 ) =-AF 
A ( 6 , 7 ) =-2 • *0M 
A ( 7 , 6 ) a2  »*0M 
TF  =  84 

NF=10.*TF  +  1.  1 
DT“ • 0 1 

DO  82  NQ= 1 ,LP 
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DO  83  LA«2»4»2 
DO  83  LB-1*7 

83  GG(LA»LB)»G3(LA.LB*NO) 

DO  62  J ■ 1 » 4 
62  X ( J  > -XX ( J ) 

DV-O. 

T»0. 

RPO-O. 

DM1-2.*(0M*T+ALF) 

X  (  5  >  -  1  • 

X ( 6 ) «COS ( DM1  ) 

X ( 7 ) "SIN ( DM1  ) 

UX-0« 

UY-O. 

DO  24  L* 1  * 7 
UX°UX-GG ( 2  tL ) *X ( L ) 

24  UY-UY-GG(4*L)*X(L) 

USQ»UX**2+UY**2 

U*SQRT(USQ)  _ 

DV»DV+U*DT 
RSQ»X( 1 ) **2+X{ 3 )**2 
R*SORT(RSQi 
RP“AMAX1 ( R  »RPO ) 

RPO^RP 

NA*  1  '  " 

V2(NA)=DV 
RP2(NA)*RPO 
R2(NA)*R 
UM ( NA ) °U 
T2(NA)=T 

WRITE  (6*14)  NQ.Q(NO) 

14  FORMAT  ( 1H1 ♦ 10X  » 3HJ  • *  1 3 » 10X * 3HQ  ■ #F8*4///1 7X * 5HPO I  NT ♦  - 
24X»4HTIME»9X*  7HDELT  A  V*6X*10HPEAK  DR  I  FT . 1 IX  * 5HDR I  F  T  * 6X ♦ 
37HCONTROL/2X )  , 

DO  54  NA*2»NF 
DO  53  NBal.lO 
CALL  RKDES  ( T *X  ♦  4 *DT *FX ) 

31  DM1=2.*(0M*T+ALF) 

X  <  5  )  - 1  • 

X(6)»COS(DMl) 

X ( 7 ) ®S I N ( DM  1  ) 

UX«0. 

UY*0. 

DO  22  L  =  1 »  7 
UX-UX-GG ( 2  » L ) *X ( L ) 

22  UY=UY-GG ( 4  *  L ) *X ( L ) 

USQ=UX**2+UY*#2 
U  =  SQR  T ( USO ) 

DV=DV+U*DT 

RSQ=X { 1 ) **2+X ( 3 ) **2 

R=SORT(RSQ) 

RP  =  AMAX 1 ( R  » RPO ) 

53  RPO=RP 
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V2(NA ) *DV 
RP2INA) »RPO 
R2( NA)  »R 
UM ( NA ) °U 
54  T 2  <  NA ) » T 

82  WRITE  (6*15)  ( L ♦ T2 ( L) * V2 ( L ) .RP2 ( L ) * R2 ( L ) *UM( L ) » L« 1 . NF ) 
15  FORMAT  (  I22.0PFQ.2, 1PE16# 7  * 1PE 16. 7. 1PE16. 7 . 1PE16. 7 ) 

S  TOP 
FND 

tIBFTC  STATE  M94.XR7 

SUBROUTINE  FX  (T»X*t)X) 

COMMON  /COM1/A ( 7.7) *  GG ( 7 . 7 ) . AT . OM » ALF 
DIMENSION  X  <  7  > »DX ( 4  ) *AX ( 4  *  7 ) 

DM1=2.*(0M#T+ALF) 

X  (  5  )  »  1  • 

X ( 6 ) aCOS (DM1) 

X ( 7 ) «S IN ( DM  1 ) 

A2 1*A( 2  *  1  ) 

A23=A( 2.3) 

A41«A( A, 1 ) 

AA3*A( 4.3 ) 

DM2=AT*X(6) 

DM3a  AT*X ( 7 ) 

A(2»1)*A(2.1) +DM2 
A(2.3)aA(2.3)-DM3 
A(A*1)=A(A»1) -DM3 
A ( 4.3) *A(A,3 1-DM2 
8  DO  3  1*1.7 

AX ( 1 »L ) *0. 

3  AX(3.L)=0. 

AX(1»2)=1. 

AX(3.A)=1. 

DO  A  Lai. 7 

AX(2.L)*A(2.L)-GG(2.L) 

A  AX(A.L)=A(A.L ) -GG  <  A . L ) 

A ( 2. 1 ) = A? 1 
A( 2.3)=A23 
A  (  A.  1  )  =  AA  .1 
4(4,3 ) =AA3 
DO  5  L=  1 , A 
DX (L ) -0, 

DO  5  K=  1  *  7 

5  DX(L)=AX(i_.K)*X(K)+DX(L) 

RETURN 

END 

SIBFTC  RUNKUT  M9A.XR7 

SUBROUTINE  RKDES ( X , Y , N » DX , F ) 

DIMENSION  Y(28), YO (28), YT(28).YP (28), P0( 28) 

DIMENSION  PI ( 28  )  ,P2 ( 28 ) .P3 ( 28  ) 

1  XO  =  X 

X  =  X  +D  X 
H=DX 

2  IF (ABS(H).GT •ABS(X-XO) )  H=X-XO 
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DO  4  I»1 »N 

4  YO< I) ■ Y (  I  ) 

HT-H 

X  T*  XO 

RMAXP  »  1.E37 
DO  5  I  *  1 »N 

5  YT ( I ) * YO ( I ) 

ASSIGN  6  TO  K 

GO  TO  20 

6  DO  7  I ■ 1 *N 

7  YP ( I ) ■ Y ( I ) 

8  HT*0.5*H 

ASSIGN  9  TO  K 

GO  TO  20 

9  DO  10  I » 1  »N 

10  YT  < I ) » Y ( I ) 

XT-XO+HT 
ASSIGN  11  TO  K 

20  CALL  F { XT » YT  *  PO ) 

DO  21  1*1. N 

21  Y(  I  )*YT<  I  )+0.5#HT#POfn 
CALL  F<XT+0.5*HT»Y»P1) 

DO  22  I - 1 *N 

22  Y( I ) -YT  < I )+.5*HT*Pl< I ) 

CALL  F<XT+0.5*HT.Y.P2) 

DO  23  1*1. N 

23  Y(I )=YT( I )+HT*P2( I) 

CALL  F ( XT+HT  »  Y  »P3  ) 

DO  24  I  °  1  .N 

24  Y(I )*YT( I ) +HT*  <  PO ( I ) +2 •* ( PI ( I  )+P2( I) >+P3( I ) )/6» 

GO  TO  K, (6.9.11) 

11  RMAX“0 • 

DO  12  I  *  1  »N 

12  RMAX *AMAX 1( RMAX.. 07#ABS( ( Y ( I )-YP( I ) )  /Y ( I )  )  ) 

IF  ( (RMAX.GT. l.E-06 )• AND. ( RMAX.LT. RMAXP) )  GO  TO  17 
XO*XO+H 

IF(XO.EQ.X)  RETURN 

I F ( (RMAX.LT.l.E-7 ) .OR. ( RMAX .GT. RMAXP  > )  H  =  H+H 
GO  TO  2 

17  H«HT 

X  T  =  XO 

DO  19  I  “  1  »N 

18  YP (I ) * YT (  I  ) 

19  YT ( I  ) = YO (  I  ) 

RMAXP  =  RMAX 
GO  TO  8 

END 

SDATA 

$N 1  AA=.7528*AB=1.267*AC=2.253»AD=.001365  . 

AE=. 002423,  AF-. 004094.  AG*. 007268.  AT=. 008393* 

0M  =  ,9252  S 

SN3  X ( 1 )  =  .002 •  X ( 2 ) =0 *  »  X ( 3  )  =0 .  .  X(4)=0.»  ALF  =  90*  . 

LM  =  4  $ 
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0*0100 

-0.72882692E  00 
0.6345B543E  00 
-0, 12758727E  01 
-0  r1724663E  00 
-0.  3760*  10E-02 
-0.50880010E-03 
-0*  5 1 252 145E-02 
0*  1  1846696E  01 
-0*  53724663E  00 
0.17579287E  01 
0.10383819E  01 
0. 19819423E-02 
-0*4277791 5E-02 
0.59098098E-02 
0.3000 

-0.48292343E-00 
0.82325358E  00 
-0.12893765E  01 
-0.46494085E-00 
-0.11963367E-02 
-0.21485960E-03 
-0.75401253E-0 2 
0.22391029E  01 
-0* 4649408 5E-00 
"  ’9210094E  01 
0*1 9406404E  01 
0*33821 66 3E-02 
-0.1 1317102E-01 
0.42906228E-02 
3.0000 

0.95098654E  00 
0.14873070E  01 
-0 • 1 199505 1 E  01 
-0. 79091 149E-01 
-0.62572880E-04 
0,27161561 E-02 
-0.96700978E-02 
0.35325503E  01 
-0* 79091 149E-01 
0.41497470E  01 
0.28243360E  01 
0.43564986E-02 
-0.12124909E-01 
-0.20369150E-02 
10.0000 

0, 25984848E  01 
0,2 1933012E  01 
-0. 14953858E  01 
0*943732 11E-01 
0.71278690E-03 
0.407336  0*1  E-02 
-0. 88293092E-02 
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0.42276491E  01 
0*94 3 732111- -01 
0.52614085E  01 
0.33002126E  01 
0.42135844E-02 
-0, 10761255E-01 
-0.32439094E-02 
tFOF 
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Performance  Calculation.  States 


$ JOR  0*15.5000  67-114, STE I NBERG » AF I T  SE 

flRJOR 

$IRFTC  MAIN  M94.XR7 

COMMON  /COM  1 /A (7*7)  *  GG  <  7  »  7  ) »AT  *OM  » ALF 
DIMENSION  H 28 ) ,X(7) ,0( 10) ,63(4,7,10) »XX(4) 
DIMENSION  V2 ( 841 ) , R  P2 ( 84 1 } ,R2(841) *T2(841)»UM(841) 
DIMENSION  X7(641),Xl(841),X2(841)*X3'841),X4(841) 
DOUBLE  PRECISION  DV 
EXTERNAL  FX 

NAMFL 1ST  /Nl/AA.AB,AC,AD,AE.AF»AG,AT,OM 
NAMEL I r  r  /N3/X,ALF,lM 
READ  (5,N1) 

READ  (5.N3) 

DO  99  J*l,4 
99  XX(J)«X(J) 

LP  =  LM+ 1 
Q ( LP ) =0 • 

DO  43  K  =  2  » 4  »  2 
DO  43  J  =  1  »  7 
4  J  G3(K,J,LP) =0. 

DO  42  NQ-l.LM 

READ  (5,13)  Q!NQ)  ,(  (G3(K»J»NQ)»Jal*7)»K.  =  2»4»2) 

42  WR I TE ( 6 , 1 3 )  Q< NQ ) , (  ( G3 ( K , J *NQ ) , J* 1 , 7  )  , K  =  2 , 4 , 2  ) 
WRITE) 6, 13)  Q(LP) ,(  < G3 < K , J , LP ) , J= 1 , 7 ) , K  =  2 , 4 , 2  ) 

13  FORMAT  (F10.4/(E17,8) ) 

ALFD= ALF 

ALF= ALF* 3. 1415927/180. 

DO  6  1=1,7 
DO  6  J=1 ,7 
6  A ( I ,J)*0. 

A  (  1 , 2  )  =  1  , 

A ( 2 , 1 ) =  AA 
A ( 2  ,3 ) =AB 
A  (  2  » 4  )  =  2  • 

A  (  2  *  5 ) = AD 
A ( 2  * 6 ) =AF 
A ( 2  *  7 ) =-AG 
A  (  3 , 4  )  =  1  « 

A  (  4, 1  )  =■  AB 

A ( 4 , 2 ) =-2. 

A ( 4 , 3 ) = AC 
A  (  4 , 5  )  ~  A  E 
A ( 4,6) =- AG 
A ( 4, 7 ) =- AF 
A ( 6 , 7 ) =-2. *OM 
A ( 7 , 6 ) =  2 • *OM 
T  F  =  84 

NF=10.*TF+1 . 1 
DT= • 0 1 
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DO  82  NO"  1  * L P 
DO  83  LA«2.4,2 
DO  83  LB-1.7 

GG(LA*LB)«G3(LA.lB,NO) 

DO  62  J=1 ,4 
X( J)-XX( J) 

DV«0. 

T«0. 

RPO=0, 

DMl«2«#(OM*T+ALF ) 

X  (  5  )  » 1  • 

X ( 6 ) “COB ( DM 1 ) 

X ( 7  )  -SIN (DM1 ) 

UX  =  0. 

UY"  0  • 

DO  24  L*l»7 
UX=UX-GG(2»L)*X(L) 

UY"UY-GG(4»L  >*X(L) 

USQ»UX**2+UY**2 
U  =  SQR  T ( USQ ) 

DV"DV+U*DT 

RSQ  =  X ( 1 ) #*2  +  X ( 3 ) **2 

R=SORT ( RSO ) 

RP= AMAX 1 (R.RPO) 

RPOaRP 
N  A"  1 

V2(NA) »DV 
RP2 ( NA ) =RPO 
R2 ( NA ) =R 
UM(NA) =U 
X7 ( NA ) =X ( 7 ) 

XI  (NA) =X ( 1 ) 

X2 (NA) =X (2 ) 

X3(NA) =X (3 ) 

X4(NA) =X (4 ) 

T2 ( NA ) =  T 

WRITE  ( 6  » 1 4  )  NQ.Q(NQ) 

FORMAT  (1H1.10X.3HJ  =,I3,iOX,3HQ  « . F 8 . 4// / 1 7X , 5HPO I  NT . 

Do’54TNA  =  j!JSF2HX7’1‘‘X'2H,<1’I‘,X’2H,<2’14)<’2HX3’14X’?H’‘‘,/JXI 

DO  53  NB=  1  ♦  10 

CALL  RKDES  ( T »X  .4 , DT  .  FX ) 

DM  1  =  2 • # ( OM*  T  +ALF ) 

X ( 5 ) =1. 

X( 6) "COS (DM1  ) 

X ( 7 ) =SIN { DM1  ) 

UX  =  0. 

UY  =  0  • 

DO  22  L  =  1 ,7 

IJX  =  IJX-GG(2.L)*X(L) 

UY  =  LJY-GG  (4  ,L  )*X(  L  ) 

USQ=UX**2+UY** 2 
U  =  SOR  T ( USO ) 
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DV«OV+W*DT 

RSQ  -  X  (  1  )  **?.  fX(  3)  **?. 

R»SQR  T ( RSQ) 

RP“AMAX 1 ( R  »  RPO ) 

53  RPO“RP 

V2 (NA) *DV 
R°2(NA) «RPO 
R2(NA)»R 
UM ( N A ) SU 
X  7  (  N  A  )  «  X  (  7  ) 

X 1 ( NA ) -  X  <  1  ) 

X2(NA)=X(2) 

X3 ( NA ) aX  (  3 ) 

X4(NA)=X(4) 

54  T2 ( N A ) s  T 

82  WRITE  (6*15)  (  L » T2 (  L ) ,X U L ) » X 1 ( L ) »X2 ( L ) * X3 < U * X4 ( L ) » 
2L*1»NF) 

15  FORMAT  ( I22»0PF8.2.  1P5E16.7) 

STOP 

END 

$ IBFTC  STATE  M94.XR7 

SUBROUTINE  FX  ( T  *  X  *  DX ) 

COMMON  /COM1/A ( 7,7)  ,  GG ( 7 , 7 ) . AT , OM , ALF 
DIMENSION  X( 7) ,DX(4 ) ,AX(4,7) 

DMl=2.*(OM*T+ALF) 

X  (  5  )  =  1  • 

X  <  6 ) “COS ( DM  1  ) 

X ( 7 ) *SI N ( DM1  ) 

A21aA ( 2  » 1 ) 

A23a A ( 2  » 3  ) 

A4 1 *A (4,1) 

A43° A (4,3) 

DM2aAT*X(6  ) 

DM3aAT*X(7) 

A(2,1)«A(2,1 ) +DM2 
A(2,3)=A(2,3)-DM3 
A ( 4 , 1 ) s A ( 4 , 1 ) -DM3 
M4»3)=A(4,3)-DM2 
8  DO  3  L= 1 , 7 

AX( 1 ,L) -0. 

3  AX ( 3 , L ) =0* 

AX( 1,2) «1* 

AX ( 3 , 4 ) = 1 • 

DO  4  L= 1 ♦ 7 

AX(2,L)=>A(2»L)  -GG  (2,1) 

4  AX(4,L)-A(4*L) -GG ( 4 , L ) 

A ( 2 , 1 ) = A2 1 

A ( 2 ,3 ) =A23 
A ( 4  *  1 ) = A4 1 
A  (  4 , 3 ) =A43 
DO  5  L= 1 ,4 
DX(  L  )  =>0. 

DO  5  K-1,7 
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RE  TURN  AX  (  L  * K  5  **  (  K  >  +r>x  (  L  ) 
END 


1C  RUNKUT  M94.XR7 
SUBROUTINE  RKDE5 ( X ♦ V »N »DX * F ) 

DIMENSION  YI28)  ,Y0(28)  »YT(28) ,YP(2b> ,P0<2«i 
DIMENSION  PI ( 28 , ,P2( 28 , f p3728  ?  °  ' 

1  XO*X 


2 

4 


5 


6 

7 

8 


9 

10 


20 

21 


22 


23 


24 

11 

12 


17 


X*X+DX 

H*DX 


IF  (  ABS  (  H  )  .GT.ABS(X-XO)  )  H«X-’XO 
DO  4  1  =  1,  N 

vom=Y(  i ) 

HT  *H 
XT  =  XO 

RMAXP  =  1.E37 
DO  5  1=1 ,N 
YTU  )=YO(I  ) 

ASSIGN  6  TO  K 
GO  TO  20 
DO  7  1=1, N 
YP (  I  ) c y (  I  ) 

HT=0. 5#H 
ASSIGN  9  TO  K 
GO  TO  20 
DO  10  1=1, N 
YT ( I  )  =Y (  1  ) 

XT=XO+HT 
ASSIGN  11  TO  K 
CALL  F(XT,YT.PO> 

DO  21  1=1, N 


Y ( I ) =YT ( I ) +0 • 5*HT*P0 ( I ) 

CALL  FI XT+O. 5*HT ,Y, Pi ) 

DO  22  1=1, N 

Y(  I ) = YT ( I ) +, 5#HT*P 1 (  I ) 

CALL  F ( XT  +  O. 5*HT »Y,P2 > 

DO  23  1=1, N 

YII) = YT ( I ) +HT*P2 ( I ) 

CALL  F(XT+HT,Y,P3 ) 

DO  24  1=1, N 

YU)=YTU)+HT*(P0(I)+2*#(PiU)+P2f 

GO  TO  K, 16, 9, 11) 

RMAX  =  0 • 

DO  12  1  =  1,  N 


I  )  )+P3(  i  )  ) /6. 


F  i  ,'qhT^  ~  '  YU  ,“YP(  I  U/YU)  )  ) 

X0  =  X0^H  06),AND#'RMAX,Lr*RMAXPM  G0  T0  17 

IF(XO.EO.X)  RETURN 

qR(|RMAX*LT* 1,E"7) *0R« IRWAX.GT. RMAXP )  )  H  =  H  +  H 


H  =  HT 


XT  =  X0 

DO  19  1=1, N 
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18  YP(  IMYT  (  1  I 
10  Y  T  (  F  )  *  Y  0  M  ) 

RMAXP  =  UMAX 

GO  ro  8 

F;  NO 

UDATA 

SN1  A A*. 7628 ,AB*1 .  26  7  ,  AO  2 . 2  5  3  , AD« 

AE*. 002423,  AF=. 004094 
OM=.9252  $ 

$N3  X (  1 ) * .004  «  X ( 2 ) a0  •  *  X ( 3 ) -0  * ,  X(4)a0»,  ALF*45., 
LM -  3  S 
0.0100 

-0 .6511705 7F  00 
0.5588241 3F  00 
-0.  1 1  179198E  01 
-0.45553420E-00 
•0. 12170604E-02 
-0.74567784F-03 
-0.42003 04 8F-02 
0.80780128E  00 
-0.45553420E-00 
0. 1251273 7 E  01 
0.68767665E  00 
0.1 3943277E-02 
-0. 19536840E-02 
0.46583379F-02 
0.0200 

-0, 70958859E  00 
0 . 60904848E  00 
-0.12282993E  01 
-0.51 3432 59E  00 
-0.13313749E-02 
-0.58913541E-03 
-0.47957293E-02 
0.10358803E  01 
-0.51 343259E  00 
0.1 5639579E  01 
0.90031803E  OO 
0.17550997E-02 
-0.32934009E-02 
0  *  55 1 74362E-0? 

0.0300 

-0.72882692E  00 
0.634585 4 3F  00 
-0.12758727F  01 
-0. 53724663F  00 
-0.  1 3  7604 10E -02 
-0. 50880010E-03 
-0.51252145E-02 
0.11846696F  01 
-0.53  724663E  00 
0.  17579287  F!  01 
0. 1038381 9F  01 
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0.19819423E-02 
-0.42777915E-02 
0  #  59098098E-02 
*EOF 


/ 
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Appendix  C 
Ana  Too  Diagrams 


The  analog  simulation  was  performed  on  two  EAI  TR  48 1 s 
coupled  together.  One  computer  alone  did  not  have  a  suf¬ 
ficient  number  of  amplifiers  or  potentiometers  to  simulate 
the  entire  system.  The  circuits  are  illustrated  In  Figures 
21  thru  24.  For  the  correct  units,  all  Indicated  variables 
are  measured  In  volts,  e.g.  4  volts  for  the  variable  lOOx 
Implies  that  x  Is  equal  to  0.04.  Time  scaling  was  un¬ 
necessary  . 


78 


ANALOG  -  BASIC  SY5TEK 


GGC/EE/67-1  0 


GGC/EE/67-  10 


81 


GGC/EE/67- 1 0 


UNCLASSIFIED 


Security  Classification 


DOCUMENT  CONTROL  DATA  •  RAD 

(ttourity  ileetlllcetlan  el  III  le,  body  e  I  ebetraet  and  indealnp  annotation  muet  tit  tnttracl  whan  the  entail  taper!  it  elf  tiled) 


I.  ORKIINATINO  ACTIVITY  (Corporate  author)  It  REPORT  ICCURI  TV  C  UAIHFICA  TION 

.  UNCLASSIFIED 

Air  Force  Institute  of  Technology  (AFIT-SE)  Ti  aaoCa - 

Wrlght-Patterson  AFO,  Ohio  45433 


}.  RtRORT  TITUS 


OPTIMAL  STATION  KEEPING  AT  THE  L4  LIB  RAT  ION  POINT 


4.  OSSCRIRTIVt  NOTH  (Type  el  report  end  Inelualve  delee ) 

AFIT  Thesis 


*.  RKPO  NT  DAT! 

June  1967 


|fl.  CONTRACT  OR  ORANT  NO. 


7#.  TOTAL  NO.  OF 

83 


7  b.  NO. 

5 


•  a.  ORISINATOR'S  REPORT  NUMPIRfJJ 


b.  PROJECT  NO. 


10.  AVAIL  ABILITY/LIMITATION  NOYICIS 


AFIT  Thesis  GGC/EE/67-10 


lb.  gi  ^mjir^rjroht  NOfJJ  (Any  ether  numbere  that  mey  be  at  at  pried 

N/A 


This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is 
unlimited. 


II.  SURRLSMSNTARY  NOTH 


II.  SRONSORINO  MILITARY  ACTIVITY 

Aerospace  Research  Laboratories 
Wrlght-Patterson  AFB,  Ohio  45433 


II.  ABSTRACT 


A  control  system  Is  devised  for  maintaining  a  space  vehicle 
In  close  proximity  to  one  of  the  earth-moon  triangular  llbratlon 
points  while  minimizing  fuel  consumption.  The  problem  Is  formu¬ 
lated  as  an  optimal  state  regulator  problem  of  variational  calculus 
and  then  modern  control  theory  Is  applied.  A  quadratic  per¬ 
formance  criterion  Is  used  which  leads  to  a  linear  feedback 
control  system.  The  feedback  gains  are  obtained  by  solving  a 
matrix  differential  equation  of  the  Ricattl  type  on  a  high-speed 
digital  computer.  Performance  of  the  resulting  optimal  system  is 
verified  and  further  analyzed  on  an  analog  computer.  (U) 


D  D  . 1 473 


UNCLASSIFIED _ 

Security  Classification 


UNCLASSIFIED 


Security  Climsiflciitlon 


Station  Keeping 
Llbratlon  Point 
Optimal  Control 
Trojan  Orbit 
Space  Technology 
Aerospace 
Automatic  Control 


_ UN 

K  A_ 

I  NOLK 

ST 

INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  Enter  the  name  and  addreia 
of  the  contractor  eubcontractor,  grantee,  Department  of  De¬ 
fame  activity  or  other  organlaatlon  (’corporate  author)  laeuing 
the  report. 

2a.  REPORT  SECURITY  CLASSIFICATION:  Enter  the  over¬ 
all  aecurity  claaelfication  of  the  report.  Indicate  whether 
"Restricted  Data"  la  Included.  Marking  la  to  be  in  accord¬ 
ance  with  appropriate  aecurity  regulations. 

2b.  GROUP:  Automatic  downgrading  le  specified  in  DoD  Di¬ 
rective  5200. 10  and  Armed  Forcea  Induatriel  Manuel.  Enter 
the  group  number.  Alao,  when  applicable,  ahow  that  optional 
marklnga  have  been  uaed  for  Group  3  and  Group  4  'as  author- 
iced. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  In  all 
capital  letters.  Titlae  In  all  caeee  should  be  unclassified. 

If  a  meaningful  title  cannot  be  selected  without  claeaifica- 
tlon,  ahow  title  classification  In  all  capitals  In  parenthesis 
Immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e.g. ,  interim,  progress,  summary,  annual,  or  final. 

Give  the  Inclusive  dates  when  a  specific  reporting  period  la 
covered. 

5.  AUTHOR(S):  Enter  the  name(a)  of  authorfs)  as  shown  on 
or  in  the  report.  Entei  last  name,  first  name,  middle  initial. 

If  military,  ahow  rank  end  branch  of  service.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  date  of  the  report  as  day, 
month,  year;  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  uae  date  of  publication. 

7 a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  l.e. ,  enter  the 
number  of  pages  containing  information. 

7b.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 
references  cited  in  the  report. 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

8b,  8c,  &  8 d.  PROJECT  NUMBER;  Enter  the  appropriate 
military  department  identification,  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9a.  ORIGINATOR'S  REPORT  NUMBER(S):  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  must 
be  unique  to  this  report. 

9b.  OTHER  REPORT  NUMBER(?):  If  the  report  haB  been 
assigned  any  other  report  numbers  (either  by  the  originator 
or  by  the  sponsor),  also  enter  this  number(s). 

10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 


imposed  by  aecurity  classification,  using  standard  statements 
such  as: 

(1)  "Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC." 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  la  not  authorised. " 

(3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
uaersaeha'  request  through 


(4)  "U.  S.  military  agencies  may  obtain  copies  of  this 

report  directly  from  DDC.  Other  qualified  users 
shell  request  through 


(5)  "All  distribution  of  thla  report  Is  controlled.  Qual¬ 
ified  DDC  usere  shall  request  through 


If  the  report  has  been  furnished  tc  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  fact  and  entar  the  price.  If  known. 

11.  SUPPLEMENTARY  NOTES:  Use  for  additional  explana- 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  ( pay¬ 
ing  (or)  the  research  and  development.  Include  address. 

13.  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  the  body  of  the  technical  re¬ 
port.  If  additional  space  is  required,  a  continuation  sheet  shall 
be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  reports 
be  unclassified.  Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  military  security  classification  of  the  in¬ 
formation  in  the  paragraph,  represented  as  (TS),  (S),  (C),  or  (U). 

There  is  no  limitation  ert  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  150  to  225  words. 

14.  KEY’  WORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report  Key  words  must  be 
selected  so  that  no  security  classification  is  required.  Identi¬ 
fiers,  such  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  an  indication  of  technical  con¬ 
text.  The  assignment  of  links,  rules,  and  weights  is  optional 


